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O 
r^ ' Abstract 

r-| ■ The method of potential envelopes is used to analyse the bound state spectrum 

-*— > ■ 

c^ ' of the Schrodinger Hamiltonian H = —A + V{r), where the Hellmann potential is 

given by V{r) = —A/r + Be~'^'^ /r , A and C are positive, and B can be positive 

or negative. We established simple formulas yielding upper and lower bounds for all 

the energy eigenvalues. 
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1. Introduction 

The Hellmann potential V{r) given by 

y (r) = -A/r + Se-^7r (1.1) 

has many applications in atomic physics and condensed-matter physics [1-11]- The 
Hellmann potential, with B positive, was suggested originally by Hellmann [1-2] and 
henceforth called the Hellmann potential if B is positive or negative. The Hellmann 
potential was used as a model for alkali hydride molecules [4]. It was used also to 
represent the electron-ion [5-6] and electron core interaction [7-8]. It has also been 
shown that the main properties of the effective two-particle interaction for charged 
particles in polar crystals may be described by this potential [9-11]. 

2. The discrete spectrum : Scaling 

Many authors have studied the eigenvalues generated by the Hellmann potential 
and have tried to estimate them [1-14]. For example Adamoski [3] used a variational 
framework to obtain accurate eigenvalues. Dutt, Mukherji and Varshni [12] and 
Kwato Njock et al [14] applied the method of large- A^ expansion to approximate 
the bound states energies. In this paper we present simple upper- and lower-bound 
formulas obtained by the use of the comparison theorem and the envelope method 
[15-18]. 

We first show that discrete eigenvalues exist for the Hellmann potential for 
all values of A > 0, S, and C > 0. This result allows us to transcend the limit 
B < A assumed to be necessary in an earlier attempt at this problem by geometrical 
methods [13]. Suppose that S < 0, then we immediately have that —{A + \B\)/r < 
V{r) < —A/r. Since both upper and lower bounds are Hydrogenic potentials with 
discrete eigenvalues, the same follows for V{r). Now we suppose that S > 0. In this 
case the concern is that, for sufficiently large B, the positive term might dominate 
the Coulomb term. We see that this does not happen by the following argument. 
The function re~'^'^ has maximum value l/(eC). Hence, for S > 0, we have 
Be~^'' jr < {B/eC)/r'^, and we conclude that -A/r < V{r) < -A/r + {B/eC)/r'^. 
But the 'effective potential' for the Hydrogenic Atom in a state of orbital angular 
momentum i is given by 

Vefiir) = -A/r + £{£ + !) /r''. (2.1) 
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Hence, again, we see that V{r) is bounded above and below by Hydrogenic poten- 
tials whose corresponding Hamiltonians have discrete eigenvalues. This establishes 
our claim. 

If we denote the eigenvalues of H = —uA + A/r + Be~'~"^ by S{u}, A, B, C) , 
and consider a scale change of the form s = r/a , and choose the scale a = uj/A, 
then it is easy to show that, 

S{uj, A, B, C) = C^uSil, 4^,^,1). (2.2) 

Hence, the full problem is now reduced to the simpler two-parameter problem 

if = -A - a/r + /3e-7r, £ = £{a,l3),a> Q. (2.3) 

3. Energy bounds by the Envelope Method 

The Comparison Theorem of quantum mechanics tells us that an ordering be- 
tween potentials implies a corresponding ordering of the eigenvalues. The 'envelope 
method' is based on this result and provides us with simple formulas for lower and 
upper bounds [16-18]. We need a solvable model which we can use as an envelope 
basis. The natural basis to use in the present context is the hydrogenic potential 

h{r) = -1/r. (3.1) 

The spectrum generated by the potential h{r) may be represented exactly by the 
semi-classical expression 

Sn£{v) = min{s + vhniis)}, (3.2) 

s>0 

where the 'kinetic potential' hn£{s) associated with the potential h{r) = —1/r is 
given, in this case, exactly by hniis) = —s^/{n + tj. If we now consider a potential, 
such as V{r) , which is a smooth transformation V{r) = g{h{r)) of h{r), then it 
follows that a useful approximation for the corresponding kinetic potential fni{s) 
is given by 

fne{s) ^ g{hni{s)) (3.3) 

If g is convex in (3.3), we get lower bounds ( ~=> ) for all n and i, and if g is 
concave we get upper bounds ( ~=< ) for all n and i . 
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For the Hellmann potential, if we use the potential h = —1/r as the envelope 
basis, then the sign of g" depends only on the sign of B. An elementary calculation 
shows that 

g"{h) = -BC'^e^^'^^/h^ = BC^rh-^"^. (3.4) 

Hence, g is convex if S > or concave if S < 0. Thus in this application of the 
method we obtain upper energy bounds for B < and lower bounds for B > 0. 
The following remarks explain briefly how these results are obtained. 

We suppose for definiteness that the transformation g(h) is smooth and convex 
i.e g" > , then each tangent to g is an afflne transformation of h of the form 

yW(r) = a(t) + 6(t)/i(r), (3.5) 

where the variables a{t) and b{t) are chosen such that the graph of the potential 
V{r) lies above the graph of the potential h{r), but it is tangential to it at a point, 
say r = t. That is to say 

V{t)=a{t)+b{t)h{t) (3.6) 

and V'{t) = h{t)h'{t). (3.7) 

This means that the 'tangential potential', V^^\r) , and its derivative agree with 
V(r) at the point of contact, r = t. 

Thus, by substituting (3.3) in (3.2), we find 

O ~ minis + g{sy^/{n + £))}, (3.8) 

r>0 

which yields an upper bound if S < and a lower bound if S > . This can be 
further simplified by changing the minimazation variable s to r by the relation, 

g{Ki{s)) = g{-s'/V{n + i)) = V{r), (3.9) 

which, in turn, implies s = (n + £)^/r^. Hence we obtain finally the following 
semi-classical eigenvalue formula involving the potential V{r) itself 

£n£^rairi{{n + if/r^ + V{r)}. (3.10) 

r>0 
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4. Results and conclusion 

We now have a simple formula (3.10) for lower and upper bounds to the eigen- 
values for the Hellmann potential. In Fig.(l) we plot the ground-state eigenvalue 
bound (full line) as a function of B for the case A = 2, C = 1, along with 
the corresponding point results of Adamowski [3] as hexagons, and some accurate 
numerical values (dashed line). It is clear from this figure that the simple approx- 
imation formula gives an accurate estimate of the eigenvlues which is an upper 
bound if S < 0, and a lower bound when S > 0, as predicted by the theory. 

If we fix A, B, and C and consider the Hamiltonian H = —A -|- vV{r), with 
eigenvalues S{v) , then from (3.10) we immediatly obtain the following explicit 
parametric equations for the corresponding energy curve {v,S{v)}, namely 

2(n + £)2 



(4.1) 
^. , {n + i)^ 2{n + i)^V{r) ^ ' 

In Fig. (2) we exhibit the corresponding graphs of the function S{v)/v'^ for B = +1 
and B = —1, again with A = 2, and C = 1, along with accurate numerical data 
shown as a dashed curve. The main point of this work is to show that by elementary 
geometric reasoning one can obtain simple semi-classical approximations for the 
eigenvalues. These results are complementary to purely numerical solutions and 
have the advantage that they are expressed analytically and allow one to explore 
the parameter space without having to attend to the arbitrary additional parameters 
and considerations which necessarily accompany numerical approaches with the aid 
of a computer. 
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Fig.(l) The eigenvalues S{B) of the Hamiltonian H = —A — 2/r + Be~^ /r for 
n = 1 and £ = 0. The continuous curve shows the bounds given by the formula 
(3.10), the dashed curve represents accurate numerical data, and the hexagons are 
the results of Adamowski [3]. It is clear that the formula provides us with upper 
bounds when B < and lower bounds when B > . 
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Fig. (2) The eigenvalue bounds (full-line) for E{y)jv^ , where £-{y) is the ground- 
state eigenvalue of the Hamiltonian H = —A + vV{r) , for A = 2, C = 1, and 
B = +1, —1, together with accurate numerical data (dashed-line) . The parametric 
equations (4.1) yield upper bounds when B < 0, and lower bounds when B > 0. 



